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ABSTRACT 

We investigate the spatial and spectral properties of non-thermal emission from clus- 
ters of galaxies at 7-ray energies. We estimate the radiation flux between 10 keV and 10 
TeV due to inverse- Compton (IC) emission, 7r°-decay and non-thermal bremsstrahlung 
(NTB) from cosmic-ray (CR) ions and electrons accelerated at cosmic shocks as well as 
secondary e ± generated in inelastic p-p collisions. We identify two main region of pro- 
duction of non-thermal radiation, namely the core (also bright in thermal X-ray) and 
the outskirts region where accretion shocks occur. We find that IC emission from shock 
accelerated CR electrons dominate the emission at the outer regions of galaxy clusters, 
provided that at least a fraction of a percent of the shock ram pressure is converted into 
CR electrons. A clear detection of this component and of its spatial distribution will 
allow us direct probing of cosmic accretion shocks. In the cluster core, 7-ray emission 
above 100 MeV is dominated by 7r°-decay mechanism. At lower energies, IC emission 
from secondary takes over. However, IC emission from shock accelerated electrons 
projected onto the cluster core will not be negligible. We emphasize the importance 
of separating the aforementioned emission components for a correct interpretation of 
the experimental data and outline a strategy for that purpose. Failure in addressing 
this issue will produce unsound estimates of the intra-cluster magnetic field strength 
and CR ion content. According to our estimate future space borne and ground based 
7-ray facilities should be able to measure the whole non-thermal spectrum both in the 
cluster core and at its outskirts. The importance of such measurements in advancing 
our understanding of non-thermal processes in the intra-cluster medium is discussed. 

Key words: acceleration of particles — galaxies: clusters: general — gamma rays: 
theory — methods: numerical — radiation mechanism: non-thermal — shock waves 



1 INTRODUCTION 

The structure and evolution of the large scale universe is a fundamental tool of cosmological investigation. Clusters of galaxies 
in particular, being the largest bound objects in the universe, have been studied extensively in order to "weight" the cosmos and 
to probe its structure. In the currently favored cold dark matter (CDM) hierarchical sc enarios struct ure formation originates 
from the collapse of primordial density perturbations seeded at an inflationary epoch JPeebleslll99"3l) . In this depiction, rich 
clusters of galaxies correspond to linear perturbations on scales of order ~ 10 Mpc comoving and, there fore, their potentia l 
wells are expected to accommodate a representative fraction of the actual mass content of the universe (IWhite et alJll993T> . 
By the same token, the environment there should be the end result of physical processes that occurred below or at that length 
scale. And, in fact, observational evidence tells us that with respect to its pristine conditions the intra-cluster medium (ICM) 
gas has been metal enriched by star formation processes, magnetized in some fashion and shock heated. 

Perhaps spurred by observational progress, in recen t years the s ubject of non-thermal process in, and their impact on, 
the ICM has attracted much attention. In this regard, IVolk et all ill 9961) pointed out that metal enrichment by galactic 
winds is likely accompanied by inter-galactic termination shocks where cosmic-ray (CR) acceleration migh t take place. In 
additi on, lEnBlin et alJ lll997l) studied the case of CR protons escaping from of radio jets into the ICM and iKronberg et alJ 
teOOll) that of bla ck hole magnetic en ergy injection into the inter-galactic medium (ICM). The issue of CR ion s produced 
at cosmic shocks jMiniati et ahllibool) and their possible dynamical role was addressed bv lMiniati et alJ i2001bJ) by means 
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of a cosmological simulation of structure formation which included direct treatment of acceleration, transport and energy 
losses of a CR component. They found that if shock acceleration takes place with some ef ficiency, CR ions, due to their long 
lifetim e against energy losses and to their efficient confinement by magnetic irregularities ijVolk et alJll99ft IBerezinskv et all 
1997 ), may accumul ate inside forming structures and store a significant fraction of the total pressure there (IMiniatil 12000: 
iMiniati et aT1l2001rl . 

In order to assess the level of non-thermal activity in groups and clusters direct observations of the associated non- 
thermal emission is required. In particular, the integrated 7-ray photon flux above 100 MeV produced in the decay of neutral 
7r-mesons generated in p-p inelastic collisions, allows direct determination of the CR ion content and the pressure support 
that they provide. For the sources of CRs mentioned in the previous paragraph, the autho rs that investigated them estimated 
a 7-ray flux above 100 MeV that should be measurable by future targeted ob servations llVolk et alJll99(i IBerezinskv et all 
ll997l:lEnfilin et al.lfl997t lBlasJfl999l: lAtovan fc Volldl20ofl Miniati et alJl2001rl . However, in a recent paper iMiniatil J2002bl) 
computed the inverse Compton (IC) 7-ray emission from shock accelerated CR electrons scattering off Cosmic Microwave 
Background (CMB) photons and for typical clusters of galaxies found it compara ble to the 7-ray emission associated to 
7r°-decay. The importance of the former proc ess was recently pointed out bv j Loeb fc Waxmanl (|200(J) i n the context o f the 
unexplained cosmic 7-ray background (see also lMiniatfc002"bllKeshet et all2003t) . According to results of lMiniatH j2O02bS), the 
7-ray flux from the le ptonic component scales with the group/cluster temperature less strongly than its hadronic counterpart 
jMiniati et al.ll2001bf) . This means that it should dominate the 7-ray emission of small structures. In addition, as we shall see 
in the ensuing sections, for conventional acceleration parameter values 7-ray from IC emission and 7r°-decay can be of the 
same order even for a Coma-like cluster. Thus, the aimed detection of 7-ray emission from clusters may not necessarily reflect 
the CR hadronic component and this should be borne in mind for a correct interpretation of the observational results. 

With the objective of a correct diagnostic of future 7-ray observational results, in this paper we compute the radiation 
spectrum and spatial distribution of non-thermal emission from groups/clusters of galaxies. We inspect the case of CR ions 
and electrons accelerated at structure shocks as well as secondary elect ron-positrons (e^). However, we neglect the CRs that 
might originate at galactic wind termination shocks llVolk et allll996l) and radio galaxies. The computed spectrum extends 
from hard X-ray (HXR) at 10 keV up to extreme 7-ray at 10 TeV. We examine the following emission processes: IC and 
non-thermal bremsstrahlung (NTB) from shock accelerated electrons, IC emission from e* and 7-rays from 7r°-decay produced 
by CR ions. We find that NTB emission is negligible throughout the spectrum. IC emission from shock accelerated electrons 
should dominate at HXR energies and be comparable to the flux from 7r°-decay at 7-ray energies. Finally, IC emission from e 
is typically below these two components at all photon energies. However, we show that all of these three emission components 
(i.e., e~, e ± and 7r°) can potentially be individually measured by future 7-ray detectors with imaging capability by taking 
into account their peculiar spatial distribution and spectral properties. In particular we find that IC emission from shock 
accelerated CR electrons dominates the emission at the outer regions of galaxy clusters, provided that at least a fraction of 
a percent of the shock ram pressure is converted into CR electrons. Therefore, a clear detection of this component and of its 
spatial distribution will allow us direct probing of cosmic accretion shocks. 

The paper is structured according to the following rationale: The numerical models for both the large scale structure and 
the CR evolution are described in &\ The results on the non-thermal emission are presented in discussed in inland finally 
summarized in 



2 THE NUMERICAL MODEL 

The objective of this study is to investigate the spatial and spectral properties of nonthermal 7-ray emission from CR electrons 
and ions accelerated at large scale structure shocks surrounding clusters/groups of galaxies. The task at hand requires accurate 
modeling of: (1) the large scale structure shocks where CR acceleration occurs; (2) the distribution within galaxy clusters 
of the baryonic gas which provides the CR targets for production of 7r° and secondary e ; and (3) the spatial propagation 
and energy losses of the CRs. Furthermore, for a consistent description of the CR ions it is important to account for (the 
acceleration due to) all the shocks that have processed the gas that ends up within the GCs themselves. That is because CR 
ions with momenta above 1 GeV/c and up to ~ 10 18 eV/c have a lifetime against energy losses longer than a Hubble time. 

This suggests that the formation of clusters and groups of galaxies be followed ab initio. Furthermore in order to reproduce 
meaningfully the shocks where the CR acceleration occurs, this should be done within the "natural" framework of cosmological 
structure formation. For this reason we have carried out a numerical simulations of structure formation which follows the 
evolution of both the gas and dark matter in the universe starting from fully cosmological initial conditions. Details on the 
cosmological simulation are given below in H2.ll In addition, with the numerical techniques described in the following sections 
( H2.21 we were also able to include the evolution of CR ions, electrons and secondary e ± by accounting explicitly for the 
effects of CR injection (at shocks and through p-p inelastic collisions), diffusive shock acceleration, energy losses and spatial 
propagation. Thus in addition to the gas distribution the simulation will also provide us with the information about both the 
spatial and momentum distribution for each CR components. In particular this information will be used in the ij^lto compute 
the spatial and spectral distribution of nonther mal 7-ray emissio n within gravitationally collapsed structures. Some of these 
results were already preliminarily presented in IMiniatil i2002alJ) . Also, the simulation upon which the study is based, has 
been used by the author to study the co ntribution of 7-ray emission from CRs in the diffuse intergalactic medium to the 
cosmic 7-ray background jMiniatill2002bl) . 

In the following subsections we describe in further detail the salient features of the employed numerical techniques. These 
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are basically the same as those used in previous related studies jMiniati et al]l2001bfl iMiniatilEiioirM . and the reader 
familiar with them can skip directly to the result section in jJ3] 

2.1 Large Scale Structure 

The formation and evolution of t he large scale str ucture is computed by means of an Eulerian, grid based Total- Variation- 
Diminishing hydro+N-body code <Rvu et allll993h . We adopt a canonical, flat ACDM cosmological model with a total mass 
density fl m = 0.3 and a vacuum energy density S1a = 1 — = 0.7. We assume a normalized Hubble constant h = H o /100 
km s^ 1 Mpc -1 = 0.67 jFreedman 2000) and a baryonic mass density, Qt = 0.04. The simulation is started at redshift z ~ 60. 
The initial density field is homogeneous except for perturbations generated as a Gaussian random field and characterized by a 
power spectrum with a spectral index n s = 1 and "cluster-normalization" as = 0.9. The initial velocity field is then computed 
accordingly through the Zel'dovich approximation. 

The choice of the computational box size is a compromise between the need of a cosmological volume large enough to 
contain a satisfactory sample of collapsed objects and numerical resolution requirements. Ideally one would like to be able 
to produce rich clusters of galaxies of the size of Coma cluster because they can be compared more easily with observations. 
However, in order for such massive objects to develop out of the initial conditions described above, the size of the computational 
box should be at least a few hundred /i _1 Mpc in size. Although this could be achieved with Adaptive Mesh Refinement 
techniques, with the available computational resources and the code employed here the above box size would impose an 
unacceptably coarse spatial resolution. Therefore the size of the computational box is set to L — 50 /i -1 Mpc. Given the 
relatively small box size the dimensions of the largest collapsed objects, with a core temperature T x ~ 2 — 4 keV, are modest. 
Therefore, when in H3.4I and H3.5I we make predictions for a Coma-size cluster (with T x ~ 8 keV) the computed simulation 
results will be rescaled appropriately. The details of the rescaling procedure will be described in those sections. 

Finally, the dark matter component is described by 256 3 particles whereas the gas component is evolved on a comoving 
grid of 512 3 zones. Thus each numerical cell measures about 100 h~ 1 kpc (comoving) and each dark matter particle corresponds 
to 2 x 10 9 r^Q. 

2.2 Cosmic-rays 

In addition to the baryonic gas and dark matter, the simulation also follows the evolution of three CR components, namely 
CR ions and electrons injected at shocks and secondary e* generated i n p-p inelastic collisio ns of the CR ions with the gas 
nuclei. This is achieved through the code COSMOCR fully described in lMiniatil (I200ll.l2002bh . 

The dynamics of each of these CR component consist of the following processes: injection, diffusive shock acceleration, 
energy losses/gains and spatial propagation. In this section we describe how each of these processes is modeled numerically 
and how it is implemented in the simulation. We notice from the outset that CRs are treated as passive quantities, meaning 
that their dynamical role is completely neglected both on the shock structure (test particle limit) and the gas dynamics. 

2.2.1 Injection and acceleration at shocks 

The first step in order to compute CR injection and acceleration at shocks, is to detect the shocks themselves. Here, as in 
previous studies, shocks are identified as converging flows (V ■ v < 0) experiencing a pressure jump AP/P above a thresh- 
old corresponding to a Mach number M=1.5. For the identified shocks both mass flux and Mach number are computed. 
This i s simply done by evalua ting the jump conditions experienced by the flow as reproduced in the numerical simulation 
(e.g., lLandau fc Lifshitzlll987l) . Shocks are found both ar ound filaments and groups/clusters. In the latter case they show a 
rather complex structure (see H8.ll and iMiniati et alJl200ft for a detailed description of cosmic shocks). An example of gas 
density distribution, velocity field and shock structure within a collapsed object is provided in Fig. which illustrates the 
dominance of asymmetry in the accretion flows and the existence of strong accretion shocks at the outskirts region followed 
by weaker ones in the inner regions. A detailed description of Fig. Q is delayed until H3.ll and we shall now focus on the 
description of the adopted CR injection/acceleration scheme. 

As we shall see in the following, the mass flux through the shock and the shock Mach number determine two important 
quantities: respectively the CR injection rate and the shape of the accelerated distribution function. Injection and diffusive 
acceleration at shocks take place on spatial scales of the order of the particles diffusion length which is 

»*>-"(<&) (ago)"' (l»r V2 - W 

where E is the particle e nergy, B i s the magneti c field strength, u s the shock speed and we have assumed for simplicity Bohm 
diffusion. As discussed in iMiniatil i200lL see also I Jones et al.l|l999Tl to properly follow injection and acceleration at shocks one 
should at least resolve not only scales ~ Ajj(p) but also the shock structure including the subshock. Because such small scales 
can not be reproduced in the current simulation, injection and diffusive acceleration are not directly simulated but simply 
modeled. In particular, given the mass flux across a shock, the assumed injection model determines the fraction of particles 
that are converted into CRs. Then, the acceleration model determines how the injected CRs should end up distributed in 
momentum space. 
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F or the injection of CR ions we follow the thermal leakage prescription dEllison fc EichleJl984l : [Questll988l : lKang fc Jonesl 
I1995T) . That is, the post-shock gas is assumed to thermalize to a quasi-Maxwellian distribution function characterized by a 
temperature T s hock- The ions/protons in the high energy tail of such distribution are assumed to be able to leak back upstream 
and undergo the acceleration mechanism. The momentum th reshold for such injec tion is set to pi n j = 2ci m p yka Tshock / rrip , 
where m p is the proton mass and c\ is a control parameter jKang fc Joneslll995D . We adopt C\ 2.5 which implies that a 
fraction ~ 10 -4 of the particles crossing the shock be injected in the acceleration mechanism. In terms of shock ram pressure, 
the pressure borne by the CR ions is a small fraction for weak shocks (M ^ 3) and reaches 30 % for moderately strong shocks 
(M ~ 4 — 6). Given our simplified injection prescription for very strong shocks [M 10) this fraction can approach unity or 
so. Therefore, for consistency in these cases the fraction of injected particles is renormalized so that the CR pressure is always 
limited below 40% of the shock ram pressure. 

As for electrons (primary e~), we simply assume that the ratio between injected CR electrons and ions at relativistic 
energies be given by a parameter R e / p . The introducti on of this paramete r simplifies the treatment of this process, which is 
quite complex and yet to be fully understood (see, e.g., lEllison et al.ll200o|. for a discussion on this issue). Observationally, for 
Gala ctic CRs this ratio appears to be in the range 1 x 10" 2 - 5 x 10~ 2 jMiiller fc Tandll987l : lMiiller fc et al]ll995tlAllen et a! " 



2001 ) . Also, based o n EGRET observational upper limits on the 7-ray flux fro m nearby clusterT l^reelamTaret^ilT^9^]R^hriei 



199S : iReimer et al]l2003f) . for a CR ion injection efficiency as assumed here iMiniatil J2002bl) constrained R e / P 2 x 10 



Therefore, in the following we set R e / P = 10 -2 and all of our results concerning primary e~ are based on this assumption. 

CR ions and electrons injected as described above, are quickly accelerated to energies much higher then the thermal 
average through the diffusive shock acceleration mechanism. As already pointed out, given the prohibitively short length 
scales over which the acceleration process takes place, its direct inclusion in the simulation is not possible. In addition, for a 
typical shock with compression ratio close to four the acceleration time scale 



W(p)=21.1 (J-) (*X 



10 2 Kms- 



yr (2) 



is much shorter than the simulation dynamical time (computational time-step). Thus, as in previous studies, we assume that 
the injected CRs are accelerated instantly and distribute them in momentum according to a power-law distribution extending 
from injection up maximum energy, i.e., 

fiP) = f0 ( ) , Pinj ^ P «S Pmax- (3) 

\Pinj J 

For CR ions /o is determined through the aforementioned thermal leackage injection model and for CR electrons it is such that 
at relativistic (^ GeV) energies their ratio to accelerated protons is ~ R e / p . According to the test particle limit adopted here, 
the log-slope of the distribution function is related to the shock Mach number, M, as q = 3(7 + 1)/[2(1 — M~ 2 )} = 4/(1 — M -2 ) 
for 7 = 5/3. This implies relatively flatter distribution functions for stronger shocks. 

Based on cluster/group properties, we found it appropriate to follow CR ions between momenta p m in = 10~ 1 GeV/c 
and pmax = 10 6 GeV/c. In fact, CR ions below p m in quick ly loose energy due to Coulomb losses and b eyond 10 r GeV/c 
their confinement within cosmic structures becomes difficult llVolk et al.lll996l : IColafrancesco fc Blasi 1998). CR electrons are 
followed between p m i n = 15 MeV and 2 x 10 TeV. The latter is a reasonable value for the maximum energy to which CR 
electrons are shock accelerated, provided a magnetic field of order 0.1/iG. 

2.2.2 Secondary Electrons and Positrons 

In addition to ions and electrons injected and accelerated at shocks, CR electrons and positrons are also produced in hadronic 
collisions of CR ions with the nuclei of the intra-cluster gas. The parent CR ions are those computed in the simulation. 
Secondary e are generated in the decay of charged muons according to 

/U — » e +v e {v e ) +v^,{v tL ) (4) 
which in turn are produced in the following reactions 

p+p ^7T ± +A, 7T ± ^^i ± +^(^) (5) 

p + p -^K ± +X, K± -> [j? + v(v) (6) 

K* -> 7T° + (7) 

Additional e* are produced in cascades analogous to these by the following interactions: p+He, a+H and a+He. These are 
included by assuming a helium number fra ction of 7.3% for the background gas and a ratio (H/ He) ~ 15 at fixed energy- 
per-nucleon for the CRs jMever et al.lll997l) . The distributions of the e thus generated are obviously related to those of the 
parent CR ions. In particular, if the latter are distributed according to a power law, so are also the former. 

2.2.3 Energy losses and spatial propagation 

The numerical treatment of the CR dynamics is complet ed by accounting for both spatial transport and energy losses/gains 
as described by the diffusion-convection equation (e.g., ISkillindll975al lbii. In order to compute a numerical solution to this 
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equation, we define a grid in momentum space by dividing it in N p logarithmically equidistant intervals (momentum bins). For 
each mesh point of the spatial g rid, Xj, the CR distribution function is then described by the following piecewise power-law 
(.Tones et aLlll999UMiniatill20Q]li 

/(x J)P ) = f 3 {^)p- q ^\ 1 < < p < Pi, (8) 

where pi... are the bins' extrema. Spatial propagation and energy losses of the accelerated CRs are then followed by solving 
numerically a diffusion-convection equation that has been multiplied by '4-7rp 2 ' and integrated over each momentum bin. 
Written in the comoving coordinates system adopted for the other simulated hydrodynamic quantities, such equation reads 
(Miniatill200lh 

9n(xj,p») 



Of 



V-[un(xj,pi)] + { + — V-u) p+fc(p) 47rp 2 /(xj,^)} P! +Q(x Jl p i ). (9) 



Here n(x,p,:) = JJ l Aitp 2 /(xj,p) dp is the comoving number density of CR in the i-th momentum bin and Q(xj,pj) is 

a comoving source term, i(xj,p), describing either shock injection or e* generated in p-p collision and integrated also over 
the i-th momentum bin. In addition, a is the cosmological expansion factor, a/a = H(z) is the Hubble parameter defining 
the cosmic expansion rate. Thus the first part of the second term in eq. describes adiabatic losses/gains. Finally, bi(p) 
represents energy losses due to Coulomb collision, bremsstrahlung, IC and synchrotron e mission for electrons and e ± ; and 
Coulomb and p-p inelastic collision for ions. All of them are fully described in Miniati (2001). In cq. © u is the velocity field of 
the baryonic gas to which the CRs are assumed to be coupled through a background magnetic field. Thus once accelerated the 
CRs are passively advected with the gas flow. Notice that eq. iPHl does not incl ude the diffusion term pre sent in the ordinary 
diffusion-convection equation. This is because, as pointed out in lMiniati i200ll see also I Jones et al]ll999l) for typical flows in 
the simulation with u ~ a few 100 km s _1 diffusion over spatial scales Ax ~ 100 kpc is much slower than advection and can 
be neglected away from shocks. This is t rue for values of the diffusion coefficients typically assumed in the literature (see, e.g., 
IVolk et aljfl99d : iBerezinskv et alJll997l) . The value inferred for these diffusion coefficients is based on the assumption that 
the magnetic field fluctuations are described by a turbulent spectrum and that the total magnetic field energy density is that 
characteristic of ^tG strong magnetic fields. 

Integrating eq. is sufficient for an accurate treatment of the ionic component. In this case, after advancing the solution 
to eq. © one time-step, both fj(x.j) and gj( x j) defining the CR distribution function in eq. JS) are computed based on the 
updated values of n(xj,p;) and by assuming that /(xj, p) is continuous at each bin interface. 

For CR electrons, however, severe energy losses and distribution cutoffs render their numerical treatment much more 
delicate. For this reason for this component in addition to n(x,Pi) we also follow the correspondent bin kinetic energy 
e(x,pi) = 4-7T p 2 f(x.,pi) T(p) dp. Here T(p) = (7 — 1) m e c 2 is the particle kinetic energy and 7 = [1 — (u/c) 2 ] _1//2 is the 

Lorentz factor. The equation describing the evolution of e(xj,pi) is obtained in analogy to eq. ©, by integrating over the i-th 
momentum bin the same diff usion-convecti on equation (more properly a kinetic equation) that has been multiplied by T(p). 
In comoving units this reads (Miniati 2 00 ll) 

9e(xj,pi) l nL , , r, , >. , 2,, ^mf^iPi ["* wn 47rcp 3 /(xj,p) 



dt 



V-[u e (x JlPi )]+ [b(p)^p 2 f( Xj ,p)T(p)] m - b{p) , P J ^ hP ldp + S{x hPi ), (10) 

J Pi _i ^m e c 2 +p 2 

where S(xj,pi) = 4ir J^' i(xj,pi) p 2 T(p) dp, b(p) now includes both adiabatic loss terms and those described by be(p), and 
the third term on the right hand side includes contributions from CR pressure work and sink terms. The evol ution of the CR 
electro ns is obtained by integrating numerically eq. @ and 11 1H through a semi- implicit scheme described in iMiniatil (200 lL 
l2002bl) . The slope of the distribution function at each grid point and momentum bin is then determined self consistently by 
the values of n(xj,pi) and e(pcj,pi). 

However, for the simulated CR electrons and e* with momenta between 10 2 GeV/c and 20 TeV cooling time due to 
IC losses, T coo i ~ 10 5 — 10 r yr, is much shorter than the typical computational time-step, At ~ 10 7 — 10 8 yr. Thus, because 
UT C ooi -C uAt ^ Ax (Courant condition), th ese particles on ly exist in grid cells where the source term i(xj,p) 7^ and will not 
propagate outside it. Thus, as discussed in <Miniadl2002bh . in this case it is computationally advantageous and numerically 
correct to take directly the steady state solution to eq. (UJ as 

/Pi PPi ^„ Z* 00 

P 2 /c(xj,p») dp = -4tt J ^-y J p 2 i(x j ,p)dp (11) 

[a similar expression holds for eq. IjlUtj. Notice that n(xj,Pi) refers to the cell- volume averaged number density of CRs. 
Physically, the expression in eq. il l 1 li represents a summation of all the individual populations of CR electrons within the grid 
zone xj that, as they emerge from the acceleration region, and are being advected away from the shock, develop a cut-off due 
to energy losses. For what follows in the next sections it is relevant to notice that for an injection spectrum i(p) oc p~ s and 
for the case of energy losses dominated by IC emission, b(p) oc p 2 , the steady state distribution given in eq. (II IP implies a 
cell- volume averaged distribution function (f(p)) oc n(xj,pi) / Ax 3 oc p~^ s+1 '. 

Finally the CR ion distribution function is mapped by 4 momentum bins. For the electrons and e" 11 we used 5+1 momentum 
bins. The first five are logarithmically equidistant and cover a momentum range between p m in = 15 MeV up to P2 = 10 2 
GeV. The last electron momentum bin stores the CR electron steady-state distribution function between p2 = 10 2 GeV and 



6 F. Miniati 



Figure 1. Color image showing the density distribution on a plane passing through the center of the collapsed object. It is in units of 
cm~ 3 and the various level correspond to values indicated by the colorbar on the right of each panel. Arrows describe the velocity field 
on the same plane. Their number has been reduced by a factor nine for clarity purposes. Finally, blue contours indicate the isolevel 
of compression (V ■ v <0) where shoks occur. Narrow countour features correspond to location of strong shocks. Left and right panel 
correspond to X-Y and Z-X coordinate planes respectively. Axis scales are measured in h~ Mpc. 

20 TeV given by eq. li lt . Our tests indicate that increasing the number of momentum bins for either CR component (ions 
and electrons) does not affect the results. 



3 RESULTS 

In the following sections we present results based on the simulation just described on the emission of high energy radiation 
from CRs in groups/clusters of galaxies. In particular, we inspect the spectra, radial dependence and synthetic images for 
various emission mechanisms and for each computed CR population. For the purpose we excise out of the simulation box 
regions containing individual virialized objects. For each of them we compute the nonthermal radiation produced through 
emission processes described below, by using the simulated hydro (density) and CR (distribution functions) quantities. For 
simplicity in the following we focus only on a couple of simulated virialized objects, although the studied radiation properties 
are independent of our peculiar choice. Finally, in H3.2I and H3.3I we describe the qualitative properties of the nonthermal 
radiation whereas in q3.4l and H3.5l we will attempt quantitative predictions for a Coma-like cluster of galaxies. 

3.1 Density and Velocity Structure 

First we consider the largest virialized object in the computational box. This is characterized by an X-ray core temperature 
T x ~ 4 keV. The object is partially described in fig. Q where the two panels show two dimensional cuts across the coordinate 
planes X-Y and Z-X respectively of the density distribution (in color scale), the velocity field (arrows) and the location of 
shocks (contours). The gas density is n c ~ 6 x 10 -3 cm -3 at the core and drops by about 2 orders of magnitudes at a 
distance of a few Mpc. The tree-dimensional velocity characterizing the accretion flows is typically of order of 10 3 km s . 
T he velocity field conv erges toward the mass concentration and becomes quite complex close to it. As described in more detail 
in lMiniati et al.l i2000h . this is due to the large asymmetry of accretion flows. Noticeably, stream along filaments, which carry 
more momentum, reach closer to the central regions before being shocked. As a result the ensuing shock structure is also 
complex. The largest shocks are found up to ~ 5 ft _1 Mpc from the core of the collapsed object. Additional shocks are present 
at intermediate distances, especially along filaments as already pointed out. These shocks have Mach numbers that range from 
up to ~ 100 in the most outlining regions, to 3 — 10 along filaments. Close to the high temperature core some discontinuities 
are also found. However, as it is apparent from the broad features of the isocontours, these are weak, tran-sonic shocks. 

For the selected object the energy in CR ions is Ecu — 6 x 10 61 erg in the inner 1.5 Mpc and corresponds to about 26% 
of the thermal pressure within the same volume. We notice the presence in the neighborhood of the collapsed objects of a 
small structure which is is marginally visible above and to the right of it, in the left and right panels of fig.Qrespectively. As 
we shall see in the following, this will affect (particularly, will increase) the computed radiation spectrum at large radii. 

3.2 Radiation Spectrum 

In fig.|5|the quantity 'e L s ' in units 'keV s -1 ' associated with the selected object is plotted as a function of photon energy, e. 
l L e ' is the volume integrated spectral power and is defined as 



where j e (e) is the spectral emissivity in units 'keV cm -3 s _1 keV -1 ' and V(R) = 4tyR 3 /3 is an integration volume defined 
by a radius R. Here we set R = 5 ft -1 Mpc in order to account for the shock accelerated electrons which are typically located 
away from the cluster center (see figQ. We consider the following emission processes: 7r°-decay (solid line), IC emission 
from shock accelerated electrons (dotted line) and from secondary e (dashed line) scattering off CMB photons, and NTB 
from shock accelerated electrons (dot-dashed line). Thermal X-ray photons, of energy ex, can be neglected as seeds for IC 
emission of 7-rays with respect to CMB photons to first order approximation. In fact, although the number of available IC 
scattering electrons is larger by a factor ((£x)/(ec«fl))' ,_3 " 2 this is not sufficient to compensate for the lower number of 
seed photons (especially for flat CR distributions such as those produced at accretion shocks) and for the IC cross-section 
suppression due to Klein-Nishina effects which enter for x ~ 1.7[e(GeV)ex (keV)] 1//2 ^ 1 (see also, e.g., lEnfilin fc Biermannl 
Il99k iBerrington fc Dermerll2002l) . 

According to fig. [5] both emission from 7r°-decay and IC from primary e~ contribute significantly to the 7-ray emission 
above 100 MeV. The actual proportion between these two components is determined by two factors: (a) The ratio of accelerated 
electrons and ions at relativistic energies, which has never been measured for cluster shocks. It is represented by the parameter 
R e / P and was set to 10 -2 in the current study, (b) The temperature of the group/cluster, from which the 7-ray luminosity 
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Figure 2. Radiation spectrum extending from 10 keV up to 10 TeV produced by the following emission processes: IC ( dotted) and NTB 
(dash-dotted) from shock accelerated CR electrons, IC emission from (dashed) and -7 -rays from ifi -decay produced by CR ions (solid). 

from 7r°-decay and IC emission depend differently. In fact, for both processes the temperature dependence is expected to be 
of pow er-law type, but with a moderately steeper index for the hadronic emission component ( Min iati et al.ll2001bl iMiniatil 
2002b). Thus the relative strength of these two contributions as presented in fig.|5]is only indicative at this stage. However, it 
is within the objectives of the present paper to investigate possible strategies for their experimental determination. At HXR 
and soft 7-ray energies below 10 MeV, IC from primary e~ dominates the spectrum although the contribution from e* is not 
negligible. However, as we shall see in the next sections, this source of emission should also be detectable at HXR energies 
owing to its different spatial distribution with respect to the primary e~ . 

Noticeably the radiation spectra of these three components are rather "flat". This indicates that the primary CR ions 
and e~'s were generated in shocks of at least moderate strength, i.e., M 4. In deed, for the hadronic componen t, the 7-ray 
emissivity from a CR ion distribution f(p) oc p~ q " , is j £ (e) tx £ -( 4 <Jp- 13 )/ 3 ( c f. iMannheim fc Schlickeiseilll993) . Thus, the 
case of an emission spectrum such that eL e oc ej e oc e°, implies a CR ion distribution with q p ~ 4. As for both primary 
electrons and secondary e* , because of the severe energy losses suffered by these particles, the emission spectrum must be 
related to their steady state distributions. As pointed out at the end of H2.2I this is characterized by a steady-state log-slope, 
q S s = qsource + 1, that is steeper by one unit with respect to the source spectrum. The produced IC emissivity is of the 
form j s (e) oc £ ~(9ss-3)/2 ^jgjj sa tisfies the condition ej £ oc e° for q a3 ~ 5 and, therefore, q SO urce ~ 4. For primary e~ this 
straightforwardly requires that these CRs were accelerated at strong shocks. For secondary e* with a source term of the form 
/ e ±(e e ±) oc e-^p- 1 ^ 3 ( c f. IMannheim fc Schlickeiseijll994h . it requires that the parent CR ions were accelerated at strong 
shocks and that their log-slope distribution be q p ~ 4. 

Finally, NTB is negligible for all purposes. Notice that latter component includes also the contribution from trans- 
relativistic electrons. These manage to propagate away from the acceleration region (i.e., a shock) for a short time before 
being re-absorbed into the thermal pool due to Coulomb losses. This trans-relativistic component enhances the NTB emissivity 
with respect to that produced by the bulk of the relativistic electrons. In fact, it dominates the emission below 10 MeV causing 
a sort of bump in the shape of the spectrum about this photon energy (dot-dashed line in fig.|2J- 

3.3 Luminosity Volume Dependence 

We now inspect the spatial distribution of the emissivity associated with the various processes presented in the previous 
section. As anticipated there, such information allows in principle a measurement of each individual emission component 
(except NTB). With this in mind, in fig.[3]we show the radial dependence of the integrated photon luminosity above 100 keV, 
that is 

L 7 (> 100 keV,R) = / L e (e,R) de (13) 

J>100 keV 

where R, the radius defining the integration volume, is now a varying parameter. Obviously, at these low energies there is no 
contribution from the CR ionic component. Thus the diagram contains only three curves illustrating IC emission from primary 
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Figure 3. Radial dependence of the integrated photon luminosity above 100 keV {or the following emission processes: IC (dotted) and 
NTB (dash-dotted) from from shock accelerated CR electrons, IC emission from e (dashed). 
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Figure 4. Radial dependence of the integrated photon luminosity above 100 MeV for the following emission processes: IC (dotted) and 
NTB (dash-dotted) from shock accelerated CR electrons, IC emission from e^ (dashed) and ifi -decay (solid). 



e~ (dotted) and secondary e* (dashed) and, for completeness, NTB from e~ (dot-dashed). This plot clearly shows how the 
radiation emitted by primary e~ and secondary e originates in spatially separate regions. In fact, the latter component 
saturates quickly within the central Mpc and dominates the total photon production there. On the contrary, the contribution 
from the former component becomes significant only at a distance of about 1.5 ft Mpc from the cluster center and keeps 
increasing up to several Mpc from there. 

The situation is analogous at 7-ray energies. This is illustrated in fig. 2] where the integrated photon luminosity above 
100 MeV, L 1 (> 100 MeV,R) - defined by an equation analogous to eq. (1131 - is plotted versus radial distance. Here it is 
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Figure 5. Synthetic map of the integrated photon flux above 100 keV (top) and 100 MeV (bottom) in units "ph cm~ 2 s~ 1 arcmin~ 2 " 
from IC emission by secondary er- (top-left), primary e~ (top-right, bottom-right), and n° -decay (bottom-left). Each panel measures 15 
h —1 Mpc on a side. 



the emission from 7r°-decay (solid line) that saturates at relatively short distances from the cluster center and dominates the 
photon emission in the inner regions. As before, IC emission from primary e~ (dotted line) reaches a substantial level only 
outside a distance ~ 1.5 /i -1 Mpc. IC emission from e* (dashed line) is now unimportant and, as in the previous case, NTB 
from primary e~ (dot-dashed line) is completely negligible. 

The radial dependence of the various radiation components presented in fig.|j|]and|I]reflects both the spatial distribution 
of the emitting particles as well as the nature of the emission process. Thus, on the one hand, both e* and tv° are produced at 
the highest rate in the densest regions where both the parent CR ions and target ICM nuclei are most numerous. Consequently, 
e* IC and 7r°-decay emissivities are strongest in the cluster inner regions and quickly fade toward its outskirts. Notice that in 
the case of fig.[3]and2]this behavior is slightly altered by the presence of a nearby object (similar to the case in fig. [SJbelow) 
which causes these integrated photon luminosities to slowly grow even after few core radii instead of leveling out. 

On the other hand, because of severe energy losses, the emitting high energy primary e~ are only found in the vicinity of 
strong shocks where they are being accelerated. Notice that in this case the shocks must be strong so that enough particles 
are accelerated to high energies to produce appreciable emission. Thus, because the strongest shocks are located a t the cluster 
outskirts rather than at its core where the ICM temperature is already high (fig. also cf. iMiniati et aljEoOol) . the spatial 
distribution of the emissivity is now reversed with respect to the previous case. Notice that, in line with this analysis, the 
drop of IC emissivity (from primary e~ ) toward the inner regions is slightly more abrupt at 7-ray energies (fig. than at 
HXR energies (fig. [3J . 



3.4 Synthetic Maps 



In this section we present synthetic maps produced by non-thermal emission processes. These maps were obtained upon 
simple integration along the line of sight of the photon emissivity in the thin plasma approximation. We use them to study 
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Figure 6. Left, Right: Synthetic map of the total photon flux above 100 keV (left) and 100 MeV (right) in units 'ph cm~ 2 s~ 1 arcmin~ 2 '. 
Left panel includes IC emission from primary e~ , secondary and NTB. Right panel includes emission due to n°-decay , IC from 
primary e~ , secondary and NTB (c). Center: Synthetic map of the bolometric X-ray emission from thermal bremsstrahlung in units 
'erg cm~ 2 s~ 1 arcmin~ 2 '. Each panel measures 15 h~ 1 Mpc on a side. 



the morphology of the emitting region. In addition we investigate which emission component contributes most to the radiation 
spectrum when the latter is extracted from different spatial regions. Unlike plots in fig. and 0] using synthetic maps will 
automatically account for projection effects. 

We will attempt to relate the present analysis to the case of Coma cluster which, for its large size and relative vicinity, 
is likely the best source candidate for the detection of non-thermal, high energy radiation. Our calculations are now based 
on a different collapsed object extracted as before from the simulation box but relatively more isolated (and therefore better 
suited for the current purpose) than the previous one. It has a temperature of a few keV, a core density of 2 x I0 -4 cm -3 . 
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This object is smaller than Coma cluster and for that reason we will renormalize the emission components as follows. 1 As 
for the hydro quantities we retain the temperature and density profiles produced in the simulation but renormalize them by 
(multiplying each of them by) the ratio of the observed and simul ated core values. The following observed values for Coma 
cluster are taken: a temperature T x — 8.2 keV IIArnaud et al.ll2001f) and a gas density of n gas — 3 x 10~ 3 cm -3 . The thermal 
energy within 1 Mpc computed according to temperature and density thus renormalized c orresponds to 1.6 x 10 63 erg. 

For the IC emission from primary e~ we have used the recent finding of lMiniatH i2002bT) indicating that the 7-ray flux due 
to this process scales with the X-ray temperature as _F 7 oc T 26 . The total number of secondary e , N e ± , and the associated 
IC flux Fj^-,, are normalized by assuming that these particles are responsible for producing Coma radio halo through emission 
of synchrotron radiation. For a given measured radio flux, S sy , and an assumed average ICM magnetic field strength, (B), 
the total number of emitting particles scales as 

Ftc^N e± cc-^, (14) 

where a ~ 1 is the radiation spectral index. The assumption that Coma radio halo is produced by secondary e* , as opposed 
to high energy e~ of different origin, allows us to also fix the total number of parent CR ions, N cr i, and the ensuing 7r°-decay 
fluxes, i ? 7r o^ 77 . In fact, since the steady-state total number of secondary CR electrons accords to 

N e ± ccN cri (p gaa ) 1 (15) 

1 + Ub/Ucmb 

based on this and eq. 1141 one finds 

1+ gfl 

F n o^ oc N crt ( Pgas ) k S sy , (16) 

where Ub = {B} 2 /8n and Uqmb are the energy density in magnetic field and CMB radiation field, r espectively. For th e 
synchrotron flux needed in 1141 and 11611 we adopted the value Siaghz = 640 mjy measured at 1.4 GHz bv lDeiss et al.Hl997l) . 
When extracting synthetic spectra in the next secti on we consider two cases for the magnetic field strength corresponding 
to (B) ~ 0.15/xG and (B) ~ 0.5/xG <Kim et al.lll990l) . The former choice implies a ratio of CR ions to thermal energy about 
30 % as opposed to almost 4% for the latter. This second case is assumed to produce the synthetic maps described below. 
Notice that its "renormalized" values would imply a lower CR ion acceleration efficiency then described in t|2.2.1l and would 
also correspond to a parameter R e / P ~ 0.05. 

It is worth pointing out that according to eq. 114H and 1II6II smaller values of (B) than assumed here would imply corre- 
spondingly higher fluxes of IC emission and 7-ray from 7r°-decay . For very weak fields, then, most of the assumed radio emitting 
particles could not be of h adronic origin due to EGRET upper limits on 7-ray flux from Coma cluster felasi fc Colafrancescd 
119991: iMiniati et al.lEoOlal . cf. also fig.®. In the opposite limit of stronger magnetic fields than assumed here, less emitting 
particles would be required to produce the same radio emission. As a consequence, the associated IC flux is reduced propor- 
tionally to the inverse of the magnetic energy density [cf. eq. 1141 1. However, as indicated by the expression in eq. 1161 . for 
q ~ 1 and (B) 2> ^/Ucmb /87t ~ 3.3 /xG, the 7-ray production from 7r°-decay reaches only a floor value. This is because even if 
the higher magnetic field implies a smaller number of radio emitting e 1 * 1 particles, in steady state a minimum rate of hadronic 
interactions is necessary to compensate for the increased synchrotron losses. In fact the enhanced synchrotron emission rate 
from the total particle distribution (oc (B) 1+a ~ (B) 2 ) is counterbalanced by the reduced number of steady state emitting e 
[oc 7V e ± ~ 1/(B) 2 - see eq. 1151 1. Thus in secondary models of radio halos, there is a lower limit to the expected 7-ray flux 
from 7r°-decay . 

In fig. we present synthetic maps of the integrated photon flux above 100 keV (top panels) and 100 MeV (bottom 
panels) for an assumed (B) ~ 0.5/iG. Right panels are associated with IC emission from primary e~ whereas left panels 
correspond to emission due to IC from secondary e* (top- left) and 7r°-decay (bottom- left). The top and bottom panels of fig. 
|S| are combined in the left and right panels of fig. HJ respectively to produce a synthetic map of the total photon flux above 
100 keV and 100 MeV. (For completeness these maps also include negligible contributions from the additional processes that 
were investigated in the previous section, e.g., NTB.) A synthetic map associated with the (not rescaled) bolometric X-ray 
emission from thermal bremsstrahlung is also shown in the central panel of fig.|S|to allow for comparison between thermal and 
non-thermal processes. Non-thermal maps are presented in units of "ph cm -2 s -1 arcmin -2 " whereas thermal X-ray surface 
brightness is in units "erg cm -2 s" 1 arcmin -2 ". 

The synthetic images in fig. show that the non-thermal emissivity is remarkably extended. In accord with the previous 
section we find that the emission from 7r°-decay (bottom left) and e (top left) is confined to the cluster core. There it creates 
a diffuse halo which rapidly fades with distance from the center. Also, as to be expected, there is a strong correlation between 
the spatial distribution of both these emission components and that of the thermal X-ray emission (fig.® - On the other hand, 
IC emission from primary e~ is distributed over a much more extended area. Moreover, it is characterized by a strikingly 
rich and irregular structure. This is a direct refl ection of the complex "web" of shocks that reside at the outskirts of galaxy 
clusters first pointed out bv lMiniati et alJ J200CT) . The morphology of the emissivity in fig. [H] and [5] is quite similar, although 
noticeably the central diffuse emission stands out more prominently in the high energy 7-ray map. 



1 Notice that this "renormalization" procedure can be applied to any other cluster for which sufficient observational data are available. 
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Figure 7. Synthetic spectra extracted from the core region (top) and the outskirts region (bottom). 



3.5 Synthetic Spectra 

In order to assess the separability of the different emission components, in fig. |7| we have produced synthetic spectra taken 
from a core (top) and an outskirts region (bottom). The extension of the core region corresponds to an angular size of 1° (or 
2 Mpc diameter at the red-shift of Coma cluster) whereas the outskirts region is defined as an annular ring with inner and 
outer radii of 0.5° and 1.5° (or 1 Mpc and 3 Mpc at the red-shift of Coma cluster) respectively. For the core region we consider 
two values for the magnetic field strength (0.15 and 0.5 ^iG). For each emission region in fig. |7]we plot the predicted flux 
contribution due to the main emission processes considered so far namely, IC emission from primary e~ (dot) and secondary 
e 1 * 1 (dash and dot-dash), and 7r°-decay (solid and long-dash). Notice that as the magnetic field strength is increased from 0.15 
to 0.5 fiG the fluxes due to IC from e 1 * 1 and 7r°-decay drop by about a factor ~ 10 (cf. ea. 1141 and lltil . Fig. Q shows that at 
high photon energy (> 100 MeV) the spectrum of radiation is indeed dominated by 7r°-decay in the core region (top panel) 
and by IC emission from primary e~ in the outer region (bottom). In the latter case, the residual 7r°-decay component is 
actually due to the presence in the field of view of a small structure north-west of the selected object (see fig. |5] since this is 
not part of the main cluster, we do not renormalize it according to various magnetic field values as it was done for the core 
emission). In principle any contribution such as this can be removed by excision of emission regions associated with thermal 
X-rays from structures appearing in the field of view. Below ~ 100 MeV the flux from the outskirts region is still strongly 
dominated by IC emission from primary e~ . The contribution from the latter is much reduced in the narrower field of view 
of the core region. However, for the larger field case it is still significant and at the level of IC emission from secondary e ± . 
The relative amount of radiation flux from these two components further depends on the actual shock structure subtended 
by the field of view. Nevertheless the point is that for magnetic field of a few tenths of /iG they are expected of comparable 
intensity. In order to separate them out one could measure the radiation flux as a function of radial distance in the outer 
regions, where it presumably is only due to CR electrons accelerated at accretion shocks, and then extrapolate its value for 
the core region. For this purpose, an angular resolution ~ 10' or so should be sufficient. 

Finally, above the spectra are characterized by a cutoff which appears at photon energies about 1 TeV. For IC emission 
from primary e~, for which this feature is sharpest, this is mostly due to the maximum momentum, p ma x = 20 TeV, of the 
accelerated particles. For the other processes, it is rather caused by absorption due to the reaction 77 — > e + e~ . In fact, pair 
production becomes important at TeV energies due to the presence of diffuse background radiation field at optical/infrared 
wavelengths. The spectra in fi g. |7| wer e thus corrected through a factor exp(— r 77 ) and by assuming an optical depth for pair 
production (e.g., ICoppi fc Aharonianll 19991) 

r 77 ( £7 )^0.14 (-ZL-) ( ^ (£52=^ fc-i (17) 

77 ^ 11 VlTeVj \2 x 10- 3 eVcm- 3 y V 0.023 ^ V ; 

wher e e« = 4mpC 4 /e 7 is the target photon energy and u(e«) is the energy density carried by the background radiation taken 
from lDwek fc Arendtl <1998ft . 
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Figure 8. Total radiation spectra extracted from the same spatial regions illustrated in fig. ^ - core for two different magnetic field 
values (dot and long-dash), outskirts (short dash) and their summation (solid and dot-dash) - compared to nominal sensitivity limits of 
future -y-ray observatories (thick-solid lines). For INTEGRAL-IBIS imagers (ISGRI and PICsIT) the curves correspond to a detection 
significance of 3a with an observing time of 10 6 s. All other sensitivity plots refer to a 5a significance. EGRET and GLAST sensitivities 
are shown for one year of all sky survey whereas for Cherenkov telescopes (MAGIC, VERITAS and HESS) for 50 hour exposure on a 
single source. 



4 DISCUSSION 

In fig.E|we plot sensitivity limits of planned 7-ray observatories together with the total radiation spectra from the core (dotted 
line) and outskirts (dashed line) regions that were presented in fig. |7| The sensitivity limits are plotted for point sources. For 
extended sources it will be worse by roughly a factor that goes as 8 source /8i n3t that is the ratio between the source size and the 
angular resolution of the instrument. For Cherenkov telescopes 9i na t ~ 0.1° In any case, according to this plot, several future 
experiments should be sensitive enough to detect the computed non-thermal emission at most photon energies. In particular 
the IBIS imager onboard INTEGRAL should readily measure the flux between 100 keV and several MeV. In addition GLAST 
and Cherenkov telescopes should be able to detect both the core 7-ray emission from 7r°-decay as well as the IC flux directly 
produced at accretion shocks above 100 MeV and 10 GeV respectively. One should notice that resolving sharp and complex 
structures such as those reported in fig. and E] is a real challenge, especially for coded mask techniques employed by several 
current and planned 7-ray imagers. Although a nominal resolution of ~ 10' in principle should be sufficient to identify some 
of the bright shock features, the final outcome of such measurements will hinge on the actual source fluxes and instruments 
performance. 

The measurement of the non-thermal radiation spectra at several photon energies spanning the range illustrated in fig. 
[7| provides important information about the physical conditions in clusters. First, besides the mere detection of HXR and 
7-ray radiation from clusters, important information is embedded in the spatial distribution of their surface brightness. In 
particular the extended emission component corresponds to the location of accretion shocks. Merger shocks have occasionally 
been observed in the core of clusters as relatively weak temperature jumps. However, strong accretion shocks have yet to 
be observed directly. Thus, detection and imaging of IC emission from primary electrons would provide an opportunity to 
directly observe these accretion shocks given their wide angular extension and their unique morphology. 

In addition, the flux about 100 keV will give us the first direct estimate of the energy density of CR electrons and, together 
with radio measurements, will allow an estimate of (or upper limits on) the ICM magnetic field. In this respect for a correct 
interpretation of the data it will be necessary to account properly for the fluxes contributed by both the CR electrons in the 
core and those directly accelerated at the outlining shocks. In fact, the former likely propagate in relatively highly magnetized 
regions and produce a substantial radio emission. On the other hand, because of the expected decline of the magnetic field 
strength toward the outer regions, the latter might generate only a weak radio emission but, nonetheless, a strong IC flux. 
Since, as shown in fig. |7| the HXR flux produced by this second component can easily dominate the total flux at this spectral 
range (the details will depend on the assumed normalizations on acceleration efficiency and average magnetic field strength) 
separating the contributions from CRs in the cluster core and external accretion shocks w ill be an important ste p for correctly 
measuring ICM magnetic fields. A similar point was qualitatively discussed already in iMiniati et al.l (|2p01a) and was also 
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addre ssed in iBrunetti et alJ (1200 J) in the context of a their "multiphase" acceleration model. In particular iBrunetti et alJ 
feOOll) pointed out that the magnetic field in the outer region, where most of the IC emission is produced due to the larger 
emitting volume, can be much lower than in the core where the radio emission originates. In that model all the radio and 
HXR emitting particles are accelerated with the same mechanism, whereas in our model the CR electrons that generate the 
HXR flux are accelerated in the outer accretion shocks and never make it to the cluster core. 

It is worth mentioning that our estimated non-thermal flux from IC emission above 20 keV is quite smaller than the 
recently reported measurements of excess of HXR radiation with respect to thermal emission. For the case of Coma cluste r 
our predictions, as illustrated in fig. E] fall short by a factor of several (iFusco-Femiano et aljll99st iRephaeli et al.lfl999h . 
An analogous estimate, ob tained for A2256 after appropriate rescaling, indicates a similar discrepa ncy by a factor ~ 30 
j jFusco-Femiano et al.l 120001) . although the upper limits on A3667 are respected by our predictions jFusco-Femiano et alJ 
l200^F ~The above discrepancies between prediction and reported detections could be readily improved by assuming that the 
electron acceleration e fficiency at shocks is larger by an order of magnitude with respect to what assumed here (see also 
lLoeb fc Waxmar]|2000l) . Interestingly, since these electrons are located at the cluster outskirts and need not be responsible 
for the radio halo emission, their hard X-ray emission would not constraint the ICM magnetic field detected through Faraday 
rotation measures. Ho wever, as it is apparent from fig. El th i s is at odds with EG RET experimental upper limits on the 7-ray 
flux above 100 MeV JSreekumar et al.lll99fJ : lReimedfl999l iReimer et al.ll2003l) . which allow for the acceleration efficiency 
adopted here to be increased by at most a factor ~ 2 IIMiniati 12002 a) . This constraint, however, holds true only as long 
as: (1) CR e~ are accelerated above momenta ~ 100 GeV/c and (2) the spectrum of the accelerated particles is not much 
steeper than computed here. These conditions, though, seem to be easily fulfilled for the case of accretion shocks. Thus, in 
this respect, if all of the HXR excess emission is diffuse in nature and not associated with individual sources, turbulent and/ 
or second order Fermi acceleration models may provide a more natural explanation because a high momentum cutoff in this 
range arises more naturally in those models (provided of course enough acceleration efficiency). 

One of the most awaited experiments is related to the measurement of the 7-ray flux at and above 100 MeV. This is 
important in order to conva l idate or rule out secondary e^ models for radio ha los in galaxy clusters (lDennisonll980t IVestrandl 
Il982l: iDolae fc Enfjlinll2000l iBlasi fc Colafrancescolll999l : iMiniati et aDl2001al) and in order to estimate the non-thermal CR 
pressure there jMlrhati^t^lTl2001br) . In this perspective the above authors estimated for near by clust e rs the 7-ray flux 
produced from the decay of neutral pions produced in p-p inelastic collisions. As pointed out in IMiniati (2002b), however, 
7-ray from IC emission can also be substantial and at the same level as that from 7r°-decay . Therefore, once again for a 
correct interpretation of the measurements, separation of these two components will be required. Estimating the 7-ray flux 
from IC emission due to shock accelerated CR electrons will be instrumental for ad dressing another issue of great interest . 
That is, the contri bution of this emission mechanism to the cosmic 7-ray background llLoeb fc Waxmanll2000l : lMiniati |[2002b; 
Kcs het et al.1120 03'1 , for which so far the only constraint is provided by an upper limit based on EGRET experiments (IMiniatil 
2002bl) . 

Photon energies of order ~ 10— 100 GeV to 1 — 10 TeV, will be investigated by Cherenkov telescopes. Such measurements 
will be complementary to that carried out at lower energies. Because the radiated energy spectrum is directly connected to 
the distribution function of the emitting particles, measuring the flux at different energies will provide information about the 
acceleration mechanism. In particular, the observed spectrum could differ from our predictions if the accretion shocks were 
modified by CR pressure. That in fact would cause t he radiation spectrum to soften at low energies and to become harde r 
toward higher energies, up to an energy cut-off (e.g., iMalkov fc Drurvll200lt feerezhko fe Ellisonlll999t iKang fe Jonedl2002T) . 
Given the different environmental conditions with respect to supernova remnants, these are interesting issues for a broad 
exploration of shock acceleration physics. Independent of the shock dynamics, measurements in this photon energy range 
should help determining the maximum momentum of the accelerated CR electrons and, perhaps, even of CR ions in case 
these maxima produce spectral cutoff before attenuation by 77 absorption becomes important. For the electronic component 
this is not excluded and could allow a clearer specification of the energy range in which they can contribute to the cosmic 
7-ray background. 



5 SUMMARY & CONCLUSIONS 

We have explored the spatial and spectral properties of non-thermal emission at 7-ray energies. For the purpose we carried out 
a simulation of structure formation including the evolution of baryonic gas, dark matter, CR ions and electrons accelerated 
at cosmic shocks as well as secondary e* generated in inelastic p-p collisions. We estimated the radiation flux between 10 keV 
and 10 TeV from CRs in collapsed structures due to IC emission, 7r°-decay and NTB and made specific predictions for the case 
of Coma cluster. We have assessed the importance of distinguishing among the contribution from different CR components 
for a correct interpretation of future experimental results and we have outlined a strategy to do so. Our conclusions are 
summarized as follows: 

• Two main regions for production of non-thermal radiation in clusters/groups of galaxies were identified: the core also 
bright in thermal X-ray and the outskirts region where accretion shocks occur. 

• The chief radiation mechanism at all 7-ray energies in the outskirts region is IC emission from shock accelerated CR 
electrons, provided that a fraction of a percent of the shock ram pressure is converted into CR electrons. A clear detection 
of this component and of its spatial distribution will allow us direct probing of cosmic a ccretion shocks. Su ch evidence would 
corroborate recent findings about extra-galactic radio emission from large scale shocks dBagchi et alJl2002Tl . 
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• In the cluster core, 7-ray emission above 100 MeV is dominated by 7r°-decay mechanism. At lower energies, IC emission 
from secondary e^ 1 takes over. However, IC emission from shock accelerated electrons projected onto the cluster core will not 
be negligible in general. 

• Separating the aforementioned emission components is important for a correct interpretation of the experimental data. 
This can be achieved in principle by measuring the spatial distribution of the detected emission. 

• Measuring the non-thermal spectrum will provide us with knowledge regarding: the injection efficiency as well as the 
ram pressure to CR pressure conversion efficiency for both electrons and ions; the energy range in which IC emission from 
CR electrons accelerated at accretion shocks can contribute to the CGB; the CR content in galaxy clusters; the dynamical 
conditions (CR modified or not) of the accretion shocks. 



ACKNOWLEDGMENTS 

This work was carried out at the Max-Planck-Institut fur Astrophysik under the auspices of the European Commission for 
the 'Physics of the Intergalactic Medium'. I am grateful to E. Churazov, T. Enfilin, M. Gilfanov, F. Aharonian and I. Susumu 
for useful discussion. The computational work was carried out at the the Rechenzentrum in Garching operated by the Institut 
fur Plasma Physics and the Max-Planck Gesellschaft. 



REFERENCES 

Allen, G. E., Petre, R., & Gotthelf, E. V. 2001, ApJ, 558, 739 

Arnaud, M., Aghanim, N., Gastaud, R., Neumann, D. M., Lumb, D., Briel, U., Altieri, B., Ghizzardi, S., Mittaz, J., Sasseen, 

T. P., & Vestrand, W. T. 2001, A&A, 365, L67 
Atoyan, A. M. & Volk, H. J. 2000, ApJ, 535, 45 

Bagchi, J., Enfilin, T., Miniati, F., Stalin, C. S., Singh, M., Raychaudhury, S., & Humeshkar, N. B. 2002, NewA, 7, 249 
Berezhko, E. G. & Ellison, D. C. 1999, ApJ, 526, 385 
Berezinsky, V. S., Blasi, P., & Ptuskin, V. S. 1997, ApJ, 487, 529 
Berrington, R. C. & Dormer, C. D. 2002, ApJ 
Blasi, P. 1999, ApJ, 525, 603 

Blasi, P. & Colafrancesco, S. 1999, Astropart. Phys. , 12, 169 

Brunctti, G., Setti, G., Feretti, L., & Giovannini, G. 2001, MNRAS, 320, 365 

Colafrancesco, S. & Blasi, P. 1998, Astropart. Phys. , 9, 227 

Coppi, P. S. & Aharonian, F. A. 1999, Astropart. Phys. , 11, 35 

Deiss, B. M., Reich, W., Lesch, H., & Wielebinski, R. 1997, A&A, 321, 55 

Dennison, B. 1980, ApJ, 239, L93 

Dolag, K. & Enfilin, T. 2000, A&A, 362, 151 

Dwek, E. & Arendt, R. G. 1998, ApJ, 508, L9 

Ellison, D. C, Berezhko, E. G., & Baring, M. G. 2000, ApJ, 540, 292 
Ellison, D. C. & Eichler, D. 1984, ApJ, 286, 691 
Enfilin, T. A. & Biermann, P. L. 1998, A&A, 330, 90 

Enfilin, T. A., Biermann, P. L., Kronberg, P. P., & Wu, X.-P. 1997, ApJ, 477, 560 
Freedman, W. L. 2000, Phys. Rep., 333, 13 

Fusco-Femiano, R., Dal Fiume, D., De Grandi, S., Feretti, L., Giovannini, G., Grandi, P., Malizia, A., Matt, G., & Molendi, 
S. 2000, ApJ, 534, L7 

Fusco-Femiano, R., Dal Fiume, D., Feretti, L., Giovannini, G., Grandi, P., Matt, G., Molendi, S., & Santangelo, A. 1999, 
ApJ, 513, L21 

Fusco-Femiano, R., Dal Fiume, D., Orlandini, M., Brunetti, G., Feretti, L., & Giovannini, G. 2001, ApJ, 552, L97 
Jones, T. W., Ryu, D., & Engel, A. 1999, ApJ, 512, 105 
Kang, H. & Jones, T. W. 1995, ApJ, 447, 994 
— . 2002, JKAS, 35, 159 

Keshet, U., Waxman, E., Loeb, A., Springel, V., & Hernquist, L. 2003, ApJ, in press, (astro-ph/0202318) 
Kim, K.-T., Kronberg, P. P., Dewdney, P. E., & Landecker, T. L. 1990, ApJ, 355, 29 
Kronberg, P. P., Dufton, Q. W., Li, H., & Colgate, S. A. 2001, ApJ, 560, 178 

Landau, L. D. & Lifshitz, E. M. 1987, Course of Theoretical Physics, Vol. 6, Fluid Mechanics, 2nd edn. (Oxford: Pergamon 
Press) 

Loeb, A. & Waxman, E. 2000, Nature, 405, 156 

Malkov, M. A. & Drury, L. O. 2001, Rep. Prog. Phys., 64, 429 

Mannheim, K. & Schlickeiser, R. 1994, A&A, 286, 983 

Meyer, J.-P., Drury, L. O., & Ellison, D. C. 1997, ApJ, 487, 182 

Miniati, F. 2000, PhD thesis, University of Minnesota 

— . 2001, Comp. Phys. Comm. , 141, 17 



16 F. Miniati 



Miniati, F. 2002a, in The Gamma-Ray Universe, ed. A. Goldwurm, D. Neumann, & J. T. T. Van, Moriond Astrophysics 

Meeting No. XXIInd, Les Arcs 
— . 2002b, MNRAS, 337, 199 

Miniati, F. 2002c, in The Universe Viewed in Gamma-Rays, ed. R. Enomoto, M. Mori, & S. Yanagita (Kashiwa: Universal 

Academy Press, Inc.), (preprint astro-ph/0212338) 
Miniati, F., Jones, T. W., Kang, H., & Ryu, D. 2001a, ApJ, 562, 233 
Miniati, F., Ryu, D., Kang, H., & Jones, T. W. 2001b, ApJ, 559, 59 

Miniati, F., Ryu, D., Kang, H., Jones, T. W., Cen, R., & Ostriker, J. 2000, ApJ, 542, 608 
Miiller, D. & et al. 1995, in Int. Cosmic Ray Conference, Vol. 3, Rome, 13 
Miiller, D. & Tang, K.-K. 1987, ApJ, 312, 183 

Peebles, P. J. E. 1993, Principles of Physical Cosmology (Princeton New Jersey: Princeton University Press) 
Quest, K. B. 1988, J. Geoph. Res. , 93, 9649 

Reimer, O. 1999, in ICRC, Vol. 4, Int. Cosmic Ray Conference, ed. . B. D. D. Kieda, M. Salamon, Salt Lake City, Utah, 89 

Reimer, O., Pohl, M., Sreekumar, P., & Mattox, J. R. 2003, ApJ, in press 

Rephaeli, Y., Gruber, D., & Blanco, P. 1999, ApJ, 511, L21 

Ryu, D., Ostriker, J. P., Kang, H., & Cen, R. 1993, ApJ, 414, 1 

Skirling, J. 1975a, MNRAS, 172, 557 

— . 1975b, MNRAS, 173, 245 

Sreekumar, P., Bertsch, D. L., Dingus, B. L., Esposito, J. A., Fichtel, C. E., Fierro, J., Hartman, R. C, Hunter, S. D., 
Kanbach, G., Kniffen, D. A., Lin, Y. C, Mayer-Hasselwander, H. A., Mattox, J. R., Michelson, P. F., von Montigny, C, 
Mukherjee, R., Nolan, P. L., Schneid, E., Thompson, D. J., & Willis, T. D. 1996, ApJ, 464, 628 

Vestrand, W. T. 1982, AJ, 87, 1266 

Volk, H. J., Aharonian, F. A., & Breitschwerdt, D. 1996, Space Sci. Rev. , 75, 279 
White, S. D. M., Navarro, J. F., Evrard, A. E., & Frenk, C. S. 1993, Nature, 366, 429 



This figure "figl.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/0303593vl 



